Demagnetized Electron Heating at Collisionless Shocks 
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Seventy measurements of electron heating at the Earth's quasi-perpendicular bow shock are an- 
alyzed in terms of Maxwellian-temperatures obtained from fits to the core electrons that separate 
thermal heating from supra-thermal acceleration. The perpendicular temperatures are both greater 
and lesser than expected for adiabatic compression. The average parallel and perpendicular heating 
is the same. These results are explained because, over the electron gyroradius, 5B/B ~ 1 and 
eScj)/Te ~ 1, so electron trajectories are more random and chaotic than adiabatic. Because density 
fluctuations are also large, trapping and wave growth in density holes may be important. 
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Collisionless plasma shocks play an important role in 
astrophysical and solar system plasmas by redistributing 
supersonic flow energy into thermal energy of plasma, 
fields, waves and accelerated particles. At the very heart 
of the physics of collisionless shocks is the question of 
dissipation, which is needed to limit the steepening of 
the shock [9j, [21 1 . A shock steady state is reached when 
the nonlinear steepening is balanced by dissipation with 
the addition of dispersion. 

The shock heating mechanisms for the ions and elec- 
trons are believed to be very different. In constrast to the 
ions the problem of the electron heating is not fully un- 
derstood, and remains controversial after some 40 years 
of research. Early theories focused on micro-instabilities 
such as ion-acoustic or lower-hybrid that should arise 
from the presence of unstable currents in the shock layer 
111 ]. The current consensus however is focusing on the 
presence of the macroscopic (DC) magnetic and elec- 
tric fields and associated potential in the so called de 
Hoffman- Teller-frame. These fields inflate and also open 
up a hole in the electron population when crossing the 
shock, if the electrons are magnetized and effectively adi- 
abatic [H, [l(l [l!|. This process is reversible and addi- 
tional scattering and/or instabilities are required to infill 
the distribution and thermalize the lower energy elec- 
trons. 

Previous studies suggested that electrons conserve fluid 
adiabaticity (T e ,±/B ~ constant, where T e> ± is the elec- 
tron temperature perpen dicular to the magnetic field) 
across the shock [Tfl, [l8|, although a statistical study 
showed some heating above that expected from conser- 
vation of the adiabatic moment Early observations 
of super-adiabatic electron heating prompted the devel- 
opment of alternative heating mechanisms [H-0] where 
electrons are non-adiabatic on a particle level and the 
electron trajectories are kinematically described. The de- 
magnetized and non-adiabatic behaviour in these models 
is due to the presence of a large-scale (macroscopic) elec- 
tric field gradient, resulting in super-adiabatic heating 
of electrons depending on the gradient scale. The shock 
transition scale size and other gradients are thus impor- 



tant not only for determining the importance of disper- 
sive effects [2 



14[, but also 



in addition to dissipation 
for determining the adiabaticity and resultant heating of 
the electrons. Non-adiabaticit y of electrons has also been 
studied in simulations were it was noted that 

different sub-populations of electrons could be over- or 
sub-adiabatic depending on injection angle. 

In this letter we use space plasma observations to show 
that the overall thermal electron heating is rarely strictly 
adiabatic but can be super- as well as sub-adiabatic (de- 
pending on Mach number), and that this is due to the 
fact that the ratio of perturbed to the average magnetic 
field, and the ratio of electric potential to thermal tem- 
perature (effectively the electric field gradient) are of the 
order of one, at all relevant length scales through the 
shock ramp transition. By introducing the concept of the 
perpendicular Maxwellian-temperature and contrasting 
it with both the conventional perpendicular moments- 
temperature and the adiabatic-temperature, we also 
show that previous studies have likely overestimated the 
amount of electron heating. We then discuss the proper- 
ties of the perpendicular Maxwellian-temperature at 70 
quasi-perpendicular bow shock crossings. 

We use data from the THEMIS multi-spacecraft mis- 
sion [l[ , during times when at least three of the THEMIS 
spacecraft were close together, and the time between in- 
dividual encounters of the Earth's bow shock was short. 

As an example, we will consider the THEMIS A 
crossing on December 7, 2010, at 14:59:26 UT (Uni- 
versal Time) near Xqse — 8.26, Ygse = —7.53 and 
^gse = 2.84 (Geocentric Solar Ecliptic coordinates), 
within about 30 seconds of similar crossings by THEMIS 
D and THEMIS E. This was a quasi-perpendicular 
(&Bn = 85°) shock crossing with super-critical Alfvenic 
and fast Mach numbers (Ma = 13, Mj? = 6). Obu is 
the angle between the upstream magnetic field and the 
normal to the shock and it was obtained from a mini- 
mum variance analysis of the magnetic field. The shock 
normal direction in GSE was n = (0.74,-0.67,-0.13). 
Multi-spacecraft timing analysis gave a shock speed of 
5.5 km/s, consistent with the above shock normal. The 
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FIG. 1: Snapshots of the electron distribution function sam- 
pled over 0.25s, for pitch angles of 75° — 105°. The snapshots 
are taken at distinctive regions across the shock. Arrows in 
Fig- E] show where the snapshots were taken, (a) Upstream 
region, (b) Behind the ramp, (c) Downstream, (d) Far down- 
stream. 



magnetic field increased from 5 nT to 30 nT across the 
shock and the plasma density increased from 30 cm-3 
to more than 300 cm-3. Quarter-second snapshots of 
typical 75° — 105° (hereafter called perpendicular) elec- 
tron distributions across this bow shock are illustrated 
in Figure 1. These distributions consist of low energy 
thermal cores and supra-thermal power law distribu- 
tions of halo electrons, originating in the upstream solar 
wind. Maxwellian-temperatures are defined as the fits 
of Maxwell-Boltzmann distributions to the thermal cores 
of the plasma and are illustrated as the red curves in 
Figure 1. These and all distributions in this paper are 
analyzed in the spacecraft frame because their transfor- 
mation to the more appropriate frame of the plasma has 
been shown to affect the Maxwellian-temperature esti- 
mate by no more than about one eV, because the plasma 
flow speed is small compared to the thermal electron 
speed. It is also noted that the spectra have been ad- 
justed for the spacecraft potential. The four snapshots 
of the perpendicular electron distributions in Figure 1 
are taken at select positions that highlight distinctive 
features of electron shock physics. The positions for 
these snapshots are indicated with arrows in Fig. [5J In 
the upstream region, panel (a), a mixture of 9 eV per- 
pendicular core electrons co-exist with halo electrons of 
higher energy and, perhaps, some shock accelerated elec- 
trons. Just downstream of the ramp in panel (b), the 
Maxwellian-temperature increased by a factor of two and 
flat-top distributions were seen. These electrons were ob- 
served with relatively little change in the thermal com- 
ponent further downstream, in panels (c) and (d). 



- 435 eV elprtrons (75 -105 ) 




FIG. 2: Shock electrons sampled by the THEMIS A space- 
craft, (a) Electron perpendicular and parallel temperatures. 
Fitting a Maxwellian to the thermal part of the electron 
distribution function gives a perpendicular (red solid) and 
parallel (blue dashed) heating of electrons that is less than 
that expected from adiabatic compression (black dashdot). 
Also plotted is the perpendicular moment-temperature that 
includes supra-thermal electrons, calculated using high reso- 
lution plasma data. Arrows indicate locations of distribution 
functions in Fig. [T] (b) Differential energy flux (DEF) for 22 
eV (thermal) perpendicular electrons (c) DEF 435 eV (non- 
thermal) perpendicular electrons. 



Time-domain data across the shock of Figure Q] arc 
illustrated in Figure [3J in which the upstream direc- 
tion is at the right of the figure and the shock ramp 
is indicated by the yellow rectangle. The black dash- 
dotted curve in panel (a) is the perpendicular adiabatic- 
temperature expected for adiabatic compression of the 
plasma. It is computed from the magnetic field data 
by making the first invariant, T e< ±/B = T ej ± >sw /B sw , 
a constant whose value was determined in the solar 
wind. The magnetic field is here sampled with 128 vec- 
tors/s by the fluxgate magnetometer. The red solid 
curve in panel (a) is the perpendicular Maxwellian- 
temperature and the green curve with diamonds is the 
perpendicular moment-temperature obtained from mo- 
ments of the spin period distribution function. The blue 
dashed curve is the parallel Maxwellian-temperature. Be- 
cause the moment-temperature includes contributions 
from the supra-thermal electrons, it is greater than the 
Maxwellian-temperature, as shown in panel (a). In the 
current data set, the moment-temperature exceeds the 
Maxwellian-temperature by factors from 1.2 to 2.0. 

Panels (b) and (c) of Figure 2 give the fluxes of 22 eV 
thermal electrons and 435 eV supra-thermal electrons. 
The thermal electrons were heated predominately in the 
shock ramp and the more energetic electrons were accel- 
erated about three seconds downstream of the ramp (as 
can also be seen in panel (c) of Figure 1). Thus, ther- 
mal heating and supra-thermal electron acceleration arise 



" ♦ ♦ ♦ * 



♦ 

♦ ** ♦ T- 



FIG. 3: Adiabicity TMaxwdiian/Tadiabatic as a function of the 
Alfven Mach number. The dashed curve is proportional to 
1/Ma- Circled points correspond to the two events presented 
in this Letter. 
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from different physical processes, so separating them, as 
is done by considering the Maxwellian-temperature in- 
stead of the moments-temperature, is necessary for a 
more careful analysis of electron heating at the bow 
shock. An important result in panel (a) of Figure 2 
is that the thermal electron heating described by the 
Maxwellian-temperature is less than half that expected 
for adiabatic electron heating in the compressed mag- 
netic field, as described by the adiabatic-temperature. 
This is the first experimental evidence that electrons may 
be heated less than adiabatically at the bow shock. We 
note that in this event also the moment-temperature is 
less than that expected from adiabatic heating. 

To determine statistically how thermal electrons 
are heated relative to adiabatic heating, 70 quasi- 
perpendicular bow shock crossings with Obu between 
46 and 89.5° have been analyzed analogously to the 
above example. The ratio of the observed Maxwellian- 
temperature to the adiabatic-temperature (called "adia- 
baticity") is plotted in Figure [3] versus the Alfven Mach 
number. It is seen that, at individual crossings, thermal 
electrons may be heated more or less than adiabatically 
by factors as great as three, although the average adia- 
baticity is 1.06 ± 0.52. The adiabaticity depends on the 
Alfven Mach number, with the adiabaticity decreasing 
with increasing Mach number. In this data set, there are 
no examples of sub-adiabatic heating for Mach numbers 
less than five and the few measurements above a Mach 
number of 10 all have adiabaticity less than one. The 
dashed curve in Figure 3 is proportional to 1/Ma and 
it describes the average behavior of the adiabaticity as 
a function of Ma, as is expected because T ei ± ~ B and 
M A ~ 1/B. 

Parallel heating has also been studied by computing 
the Maxwellian-temperatures for electrons having pitch 



FIG. 4: Dimensionless perturbation ratios of magnetic field, 
electric field and density for time scales r around the con- 
verted electron gyrodiameter. Quantities are measured with 
128 samples/s. (a) Total magnetic field, (b-c) 8B/B. (d-e) 
8E ■ 2p e /T e ,j_ ~ 5E/2Q. (f-g) 8n/n. 



angles between 0° - 30° and 150° - 180°. The average 
[(parallel heating) /(perpendicular heating)] is 0.95±0.14, 
which means that the parallel and perpendicular elec- 
trons are generally heated by the same amount when 
crossing the bow shock. 

The conservation of the first adiabatic invariant, 
T e _±/B, is derived under the assumption that the mag- 
netic field does not change over the gyration orbit of 
an electron or the time for it to complete one gyra- 
tion. Also assumed is that the electric potential across 
the electron gyrodiameter is small compared to the elec- 
tron thermal energy, i.e. that the electric field gradi- 
ent is small. To test the validity of these conditions, 
we compute the fluctuations of the magnetic field, B, 
and electric field, E. In this case the magnetic field re- 
quirement for conservation of the first invariant becomes 
5B/B(t, t) = [B{t + t) - B(t)]/[B(t + r) + B(t)}/2 < 1 
where r is a time scale, here corresponding to the con- 
verted gyrodiameter. 

In Figure |4j 5B/B is calculated for a randomly se- 
lected bow shock crossing by THEMIS E, for which SB 
was determined from search coil data and B came from 
fluxgate data. The local time of the crossing was 12:40 
and the latitude was 11°, with the characteristic shock 
parameters Ma — 8.8, 0g„ = 64°, and plasma (3 = 3.5. 
Panel (a) of this figure gives the magnetic field strength 
and panels (b-c) give the time domain values of SB/ B for 
the time scales, r, associated with each panel. For T e< ± 
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= 20 eV, B = 20 nT, and a shock velocity of 5 km/sec, 
the electron gyrodiameter is about 1 km and the time for 
the shock to move one gyrodiameter is about 0.2 seconds. 
Over a typical shock ramp the gyrodiameter may change 
by a factor of about two, so appropriate time scales, cor- 
responding to the shock traveling one electron gyrodiam- 
eter, are chosen accordingly. For these time scales there 
are substantial intervals in the ramp where SB/B is the 
order of or larger than unity. The SB/B is on average 
large for scales ranging from c/uj pe to pi (not shown). 
These scales are spatial in the electron frame, when the 
thermal speed of these electrons is much greater than 
the phase speeds of whistler or lower-hybrid waves in 
this frequency regime. Finite amplitude fluctuations on 
shorter timescales (temporal) will also break adiabatic- 
ity. Thus, electrons cannot obey adiabatic conservation 
at such crossings. 

A similar calculation may be done for the electric field 
fluctuations. Here the requirement is that the voltage 
across the electron gyrodiameter be small compared to 
the electron temperature, or SE ■ 2p e /T ej _i_ <C 1 for the 
electrons to be adiabatic, where p e is the electron gyro- 
radius. For T e ,j_ ~ 20 eV and 2p e ~ 1 km, T e ^±/2p e is 
20 mV/m. So this requirement becomes the condition 
that 8E/2Q be small for the electrons to have adiabatic 
trajectories. Figure HI panels (d-e), show 8E/2Q for the 
same time scales as for the magnetic field. The electric 
field component in this calculation is the normal compo- 
nent found from a minimum variance analysis. The elec- 
tric field fluctuations are comparable to or greater than 
one, so the electrons cannot move adiabatically in such 
a fluctuating electric field. It is interesting to note that 
the magnetic field fluctuations of Figure [4] dominate near 
the foot of the ramp and the electric field fluctuations are 
most important near the top of the ramp. These fluctua- 
tions are not bipolar electrostatic structures Q at debye 
scales in the ramp. We also observe such structures but 
they are intermittently spaced, while the electron heat- 
ing happens continuously throughout the ramp, which 
suggest that the electrons are scattered in the turbulent 
fields. We also note that the density fluctuations in Fig. [4] 
(f-g) are large, Sn/n ~ 1, so particle trapping and wave 
growth in density holes could be important. 

In conclusion, due to the large amplitude fluctuations 
of both B and E, the electron trajectories must be more 
random and chaotic than orderly and adiabatic. The 
thermal heating of electrons is therefore likely to be of 
stochastic nature. This result explains both the non- 
adiabatic heating of the perpendicular electrons and why 
the parallel and perpendicular heating are comparable. 
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